Simulating quantum effects of cosmological expansion using a static ion trap 
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We propose a new experimental testbed that uses ions in the collective ground state of 

a static trap for studying the analog of quantum- field effects in cosmological spacetimes, 

including the Gibbons-Hawking effect for a single detector in de Sitter spacetime, as well 

f— ^ as the possibility of modeling inflationary structure formation and the entanglement signa- 

,-H ture of de Sitter spacetime. To date, proposals for using trapped ions in analog gravity 

^^ experiments have simulated the effect of gravity on the field modes by directly manipulating 

f^ the ions' motion. In contrast, by associating laboratory time with conformal time in the 

toX) simulated universe, we can encode the full effect of curvature in the modulation of the laser 

^ used to couple the ions' vibrational motion and electronic states. This model simplifies the 

■^^ experimental requirements for modeling the analog of an expanding universe using trapped 

-vj ions and enlarges the validity of the ion-trap analogy to a wide range of interesting cases. 

I 1 PACS numbers: 03.65.-w, 04.62. -f-v, 37.10. Ty, 42.50.Ct 

4i I. INTRODUCTION 

^ , Curved-spacetime quantum field theory (QFT) is an essential ingredient of our modern under- 

I standing of cosmology [Ij and has also given rise to remarkable insights and puzzles concerning 

the nature of Planck-scale physics [2] . Unfortunately, most predictions of curved-spacetime QFT 

►j^ are extremely difficult to test directly, confined as we are to a low-curvature region of spacetime, 

■^ while dimensional analysis conspires to make the most straightforward effects extremely tiny. This 

£_J situation has given rise to a number of proposals for testing analog curved-spacetime QFT effects, 

(^ where spacetime is replaced by a system of coupled atoms, and field modes are replaced by the 

ly-T collective normal modes of oscillation of the atoms — effectively a "phonon field." The advantage 

^^ is that while the observable effects are conceptually the same (analogous), the parameters of the 

^^ experiment may be adjusted to result in a much stronger and more easily observed effect. The 

. . most common analog curved-spacetime QFT proposals involve BECs, liquid helium, or ion traps 

. 5^ as the analog spacetime O |3| . 

^ Here, we employ a general but often overlooked concept to extend the range of possible analog 

spacetime QFT experiments in ion traps, and we discuss several analog cosmology experiments 
that appear to be greatly simplified by the new picture, proposing one of them (simulation of 
the Gibbons-Hawking effect [5]) in detail. The main idea, which has been proposed previously 
in the context of BEG simulations f6| [3, is to encode the effects of spacetime curvature into 
a modulation of the analog detectors, rather than in the motion of the atoms which comprise 
the analog spacetime. We will show how this amounts to making a shift in the analogy away 
from laboratory time -f-)- detector proper time and towards laboratory time O conformal time. In 
general, such a shift in analogy can be made between the experimental lab time and any new 
time coordinate in the simulated spacetime, but in the case of Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmologies in the conformal vacuum, we will see that this particular choice 
(laboratory time made analogous to conformal time) encodes the entire effect of the spacetime 
curvature and therefore enables us to propose analog curved-spacetime QFT experiments in which 
the analog spacetime is fixed (resembling fiat Minkowski space), thus making our proposals far 
more accessible experimentally. 



^ 



Building on efforts to implement elementary quantum gates, ion trapping now offers some of the 
most exquisite demonstrations of quantum control in modern physics. The ability to couple the 
internal electronic state of one or more trapped ions to the collective vibrational degrees of freedom 
enable the ability to produce a large class of pure quantum states [HHTU]. including Fock states, 
squeezed states, and cat states, as well as quite complex information processing tasks [TT|, [12]. 
We use the control of the trapping potential afforded by ion traps, together with the ability to 
reach quantum limited motion, to propose the construction of new ion-trap analogs of the single- 
detector response in FLRW cosmology (including the Gibbons-Hawking effect [5j), along with 
several proposals for ways to apply the model to situations involving multiple detectors, including 
inflationary de Sitter cosmology |13tll4j and the entanglement signature of de Sitter spacetime |15] . 
We discuss the most important experimental parameters governing the simulation and the regime in 
which the analog experiment is likely to conform closely to the simulated cosmology. Highly efficient 
quantum measurements, based on fluorescent shelving [5], provide a practical means to test our 
predictions. The main advantage of trapped-ion analog gravity experiments over other proposals, 
such as BECs |6l [Tj, lies in the exquisite quantum control over state preparation, interaction, 
and readout. In particular, the fact that each ion is equipped with its own highly efficient and 
tunable motional detector (through laser coupling [TO]), which counts individual phonons, means 
that ion traps have a particular advantage in efficient quantum control of the detector interaction 
and readout [3j. We compare our proposal to previous work [16, 17J that suggests using ion traps 
with a time-dependent trapping potential to simulate cosmological particle creation. 

We organize the paper in the following way: Section [H] describes in more detail the conceptual 
shift from the standard "field picture" of analog QFT to the "detector picture," which we use 
extensively, and reviews the relevant properties of FLRW cosmology that make it possible. In 
Section |III[ we review the fundamental features of ion-trap simulations of cosmology, which are 



essential for describing our proposals. Section IV is devoted to the description of our experimental 



proposals. Our concluding remarks are contained in Section |Vj Unless otherwise specified, we use 
natural units, in which h = c = 1. 

II. CONFORMAL INVARIANCE: FIELD PICTURE AND DETECTOR PICTURE 

To build an experimental analogy to quantum field effects in spacetime, we will make use of two 
basic ingredients: (1) a coupled, discrete system of oscillators, where deviation of a given oscillator 
from its equilibrium position Qnit) plays the role of the field (/)(x, t), and (2) one or more detectors, 
which can be coupled to the analog field in a way which models a given type of particle detector in 
spacetime. In the case of the ion traps we discuss in the next section, for example, the electronic 
state of a specific atom is coupled to the deviation from its equilibrium position. 

Before we discuss the analog experiment, however, we recall a few basic properties of standard 
Unruh-DeWitt detectors in FLRW spacetime [18] . Recall that the usual coupling of the field to 
detector is expressed via the following interaction Hamiltonian (which, importantly, assumes that 
the cosmic time variable t equals the proper time r for the detector): 

Hintit) = f{t)m{t)^{^,t) , (1) 

where all operators are in the interaction picture, m{t) is the detector monopole moment, and 
f{t) is a "window function" modeling a time-dependent coupling of the field to the detector. In 
the case of a single detector with energy gap E at a fixed location x, this leads to a detector 
response function given by 

/OO /'OO 

dt dt'f{t)f{t')e'^^'-''W+{t,t'). (2) 

-OO J — OO 



This is the probabiUty that the detector, starting in its ground state \g), will become excited 
into its higher-energy state |e) after interacting with the field through application of the Hamil- 
tonian Hi^tit) over a very long time — i.e., at least long enough to cover the (nontrivial) support 
of f{t). This formula assumes the interaction is weak enough for time-dependent perturbation 
theory to be used. Here, D~^{t,t') is the two-point function given by 

^+(t,O = (</'(x,t)0(x,t')) (3) 

for a given state of the field (recall that the detector is fixed at spatial position x), and the bar is 
used with foresight to indicate that the time coordinate being used is t [Cf. Eq. ([9]), below]. 

An analog experiment is one which reproduces the essential features of the response function in 
curved spacetime for an analogous detector coupling and state of the oscillators (typically with the 
ground state of the collective normal modes of the oscillators being analogous to a vacuum state of 
the spacetime field). What we would like to emphasize here is that this can be achieved in multiple 
ways: by manipulating the collection of oscillators themselves or by manipulating the detector. 
When the spacetime involved is conformal to Minkowski spacetime, as in FLRW cosmology, the 
entire effect of curvature can be encoded completely within either the oscillators or the detector, 
leading to two distinct "pictures" for analog experiments. 



A. The Field Picture 

In the field picture, the relevant features of the analog spacetime are encoded in the os- 
cillators and thus directly in the modes of the analog field and the resulting two-point func- 
tion. Thus, in the response function, the analogy is contained in the two-point function di- 
rectly, i.e., ((/>(x, t)(/)(x, f')) •(-)• {4)n{t)(t)n{t')) (where 4>n{t) is an appropriately dimensionless version 
of quit))- The analogy is completed by physically moving the oscillators — an expanding array of 
ions simulating an expanding spacetime, for example |16|ll7j. The important feature of this picture 
is that t labels both the proper time of the analog detector in the lab and the proper time of the 
detector being simulated. 

The field picture is thus immediate conceptually, and for this reason it has been the "default" 
picture for proposed analog cosmological QFT experiments in the literature. We emphasize, how- 
ever, that it is not the only picture for cosmological models and that adequate control of the analog 
spacetime (the oscillators) often presents a set of significant experimental challenges. 

B. The Detector Picture 

When a given spacetime is conformal to Minkowski spacetime, and when the fields considered 
are massless and conformally coupled, we can switch to the detector picture, in which the array 
of oscillators is held fixed (in analogy to Minkowski spacetime), but the detector is modified in 
a time-dependent fashion to create the same analog effect as a curved spacetime. The important 
conceptual feature of this picture is that the laboratory proper time describing the evolution of 
the analog detector is not analogous to the proper time of the simulated spacetime detector. 
Instead, the experimental proper time (i.e., the time measured by the clock on the wall in the 
laboratory) is analogous to the conformal time of the simulated spacetime and detector. This is 
accomplished by introducing a time-dependent energy gap in the analog detector, which simulates 
an ordinary, constant-energy gap of a detector in a curved cosmological spacetime. The payoff^ for 
the conceptual shift and the modulation of the detector is a greatly simplified treatment of the 
analog spacetime — the array of oscillators is held fixed. 



To see how this works, we recall the transformation properties of fields in spacetime (see Ref. 
for more details). For simplicity, we work with a scalar field. Recall that a spacetime (with 
metric 5^^) conformal to Minkowski spacetime (wih metric hfj^^) can be coordinatized such that 
9'/ii/(x, C) = ^(x, Q'^h^y. (In this section x = (x, Q), with Q playing the role of a time coordinate.) 
In these coordinates, the two-point function in the conformal vacuum satisfies 

Z)+[x(C),C;x'(C'),C'] = ^(x)('-")/'i?^[x(c),C;x'(C'),C']f^(^')^'""^/', (4) 

where D^[x.{C,),C,]->d{C,'),C,'] is the standard Minkowski two-point function for a scalar field along 
the trajectory x((^), n is the number of spacetime dimensions, and we omit the bar over D^ because 
these are functions of C, instead of t [Cf. Eq. ^ above and Eq. ^ below] . The variable C, is thus 
interpreted as the conformal time of an FLRW universe on the left-hand side of the above equation 
and the Minkowski time on the right-hand side. 

In FLRW cosmology with no spatial curvature, the metric in comoving coordinates takes the 
form 

n-l 

ds'^ = df-a{tf"^{dx'f, (5) 

where t is cosmic time and also the proper time for a comoving observer (i.e., an observer on an 
inertial trajectory of constant x). Transforming to conformal coordinates by substitution of the 
conformal time C,, which satisfies 

dC = ^y (6) 

a{t) 



results in a metric of the following form: 
where 



n-l 

i\2 



de-Y,{d 



X 
i=l 



(7) 



^(0 = a[t(C)] . (8) 

By using Eq. Q, the two-point function on the trajectory of a comoving detector at a fixed 
location x can be expressed as 

^FLRW(*^ t') = ^FLRW(C, C) = J^(C)('-")/'^(C')^''"^/'^lIl(C, C) , (9) 

where D+{C,C) is short for D+[x{C),C;x' {('),('] for both the FLRW and Minkowski Green's 
functions as functions of C, and we explicitly show the connection to the Green's function D'^{t, t') 
used in Eq. Q. The factor of r2(C)(^~")/^f](C')(^~"')/^ may be absorbed into the definition of the 
detector window functions, and thus an analog Minkowski two-point function is equally good for 
describing FLRW cosmology, provided we can re-interpret the experimenter's "clock time" as the 
conformal time (. For this interpretation to work, however, the detector must be designed to evolve 
with respect to the simulated proper time t as e~*^* in order to model a comoving detector, rather 
than evolving in the simulated conformal time C as e"*^''. 

The appropriate modification of the detector requires that we simulate time evolution in t 
with a t-independent Hamiltonian Ho by means of a (-dependent Hamiltonian Ho{C)- This time- 
dependent Hamiltonian will be the physical Hamiltonian acting on an individual oscillator in the 



lab [HI [71 [19]. To find Hq{(), we simply write the desired Schrodinger equation in t and transform 
to C- Thus, a time-independent Hq in the equation 



i^* = Ho^ (10) 



can be re-expressed as 



i^^ = Ho{C)^, (11) 

where -ffo(C) — w^o- Simulating a detector with fixed energy gap E in t, given by Hq, requires 
a time-dependent "laboratory Hamiltonian" Hq{(^), which corresponds to a simple time-dependent 
scaling of the energy gap in laboratory time. For the case of a FLRW cosmology simulation, this 
is given by 

HoiO = ^{OHo = a[t{C)]Ho . (12) 

Observe how these ingredients come together to describe the response of a comoving detector 
at a fixed location x in an FLRW universe: 

/OO /"OO 

dt dt'f{t)f{t')e'''('-''W+^^^{t,t') 

-OO J —OO 

C(i=oo) p((t=oo) ij. jj./ 

C(t=-oo) Jcit=~oo) «(, dC, 

C(t=oo) K{t=oo) 

dC / <F(C)F(C')e^^W^)-*'(^')lD+ (C, C) , (13) 

({t=_oo) J(;{t=-oo) 



where 



F{C) = ^^(C)^''")/V[i(C)] = a[t(C)]('-")/V[t(C)] (14) 

is the window function in (^ (conformal time, also Minkowski time) necessary to simulate the window 
function f(t) in cosmic time. Note the non-standard exponential factor in the last line, coming 
from the C-dependent energy gap. This connection between the field picture and the detector 
picture is illustrated schematically in Figure [l} 

To summarize, the behavior of a detector with fixed energy gap E, determined by Hq, on a 
comoving trajectory in an FLRW universe in the the conformal vacuum may be simulated by a de- 
tector with a time-dependent energy gap, determined by Hq{C,), moving inertially in the Minkowski 
vacuum. In addition, the detector's window function must also be modified: f{t) i— >• -F(C)- The 



relevant equations for these substitutions are Eqs. (12) and (14), respectively. This correspondence 
applies equally well to multi-detector effects, such as entanglement, and it may be exploited in the 
contruction of analog gravity experiments. 

III. COSMOLOGY IN AN ION TRAP 

With the power of the detector picture, a vast number of gravitational models become ex- 
perimentally accessible simply by modulating the detector parameters appropriately. Our analog 
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(b) Detector picture 



FIG. 1. Schematic representation of the field picture and detector picture for a detector in an expanding 
universe exhibiting conformal invariance. In the field picture (a) , expansion of the universe (shown left to 
right) causes the field modes to redshift, while the detector resonant frequency remains unchanged. The 
field modes are defined with respect to cosmic time t, using the metric from Eq. (pi), and the detector Hamil- 
tonian Hg is independent of t, satisfying Eq. (10 1. Conformal invariance allows an equivalent description in 



the detector picture (b) , in which expansion of the universe (shown left to right) causes the detector resonant 
frequency to blueshift, while the field modes remain unchanged. The field modes in this picture are defined 
with respect to conformal time C, using the metric from Eq. (l7|, and the detector Hamiltonian Ho{(), given 
by Eq. (Il2]), is now C-dependent, satisfying Eq. (111. With additional modulation of the detector window 
function [Eq. (14|], both pictures result in exactly the same prediction for the physical response of the 
detector, as shown in Eq. (13). 



gravitational model of choice is the linear ion trap ^lOj . Previous authors have considered ion traps 
for simulating de Sitter spacetime [16] (including the Gibbons-Hawking effect [5]) and cosmological 
particle creation due to a finite-duration period of expansion |17| . We will compare our proposal to 
these in due course. We also note that trapped ions rotating around a ring trap have been proposed 
for detecting the acoustic analog of Hawking radiation [20j . In what follows, we describe the ion- 
trap analogs of the field (in the Minkowski vacuum) , an ordinary De Witt monopole detector |21] , 
and a modulated detector. 

We will go into some detail in describing the analogy between trapped ions and fields in space- 
time in order to make certain nuances clear, which can have an important effect when generalizing 
the model to backgrounds other than Minkowski spacetime. The reader is encouraged to notice, 
specifically, that the actual analog of the field (/> is a dimensionless discrete field </>„ [Eq. (19)], 



that the usual way of writing the ion-laser interaction [Eq. (|22[)] conflates interaction-picture time 



dependece [Eq. (20)] with additional time dependence from the laser, which has important conse- 



quences when trying to apply Eq. (12), that the Lamb-Dicke parameter rj [Eq. (23)] becomes time 
dependent when the trap frequency changeajand when the laser frequency changes [Eq. (29)], and 
that changing the laser frequency may, in some cases, also affect the Rabi frequency Qq, again 



through Eq. (29). We restore explicit use of h in this section. 



A. Discrete phonon field 



In a linear ion trap, a string of up to about 10 ions is confined to one dimension using a combi- 
nation of radio- frequency and static electric fields and capped at the ends by static electrodes [8]. 
Linearizing the total potential created by the harmonic potential due to the trap electrodes and 
the mutual Coulomb repulsion between the ions, we can model the collective motion of the ions in 
a linear trap as a collection of coupled harmonic oscillators. With a coupling matrix A obtained 



^ While this does not affect our proposal, it does affect Refs. |16II17) and any other field-picture proposals for ion-trap 
simulations of gravitational effects. 



from the linearized combination of these potentials, the Hamiltonian is given by |10j 

where q = (gi, . . . , q^) i^ ^ column vector of operators corresponding to the displacement of each 
ion from its equilibrium position, p = (pi, . . . ,Pn)'^ is a column vector of corresponding momenta, 
M is the mass of each ion, and z^ is the effective harmonic trap frequency provided by the trap 
electrodes (typically, z^ ~ a few MHz [8]). The coupling matrix A is symmetric and positive- 
definite. Thus, we can diagonalize it as A = B AB, where A = diag(^i, . . . ,//Ar) is a diagonal 
matrix of positive eigenvalues, arranged in ascending order. The orthogonal matrix B defines the 



transformation to normal-mode coordinates Q = Bq and momenta P = Bp. Equation (15) can 
be rewritten in these coordinates as 

^ = E 5f + ^^p = E H 4«p ' (16) 



2M _ 



1^^ (17) 



where 



is the oscillation frequency of normal mode p, and (since the normal modes oscillate independently) 
we have diagonalized this Hamiltonian using the standard prescription for a set of independent har- 
monic oscillators, ignoring the zero-point energy. The normal- mode raising and lowering operators 



satisfy [ap,al,/] = dpp/, as is appropriate for independent bosonic modes. The local oscillations q. 
about the ions' equilibrium positions are given, in the interaction picture, by pU 



m 




where m S {1, . . . , A^} labels the ion, and bm is an entry of B, which defines the spatial mode 
functions of the normal modejj We can define a unitless version of this local displacement oeprator, 
as well: 

N Ap) 

4>m{0 ■■= Y. -^S^'^^'^'-^V + e^'^''^^) ' (19) 

p=l ^P 

with qm{C) = ^/h/2Mi'(l)m (C) ■ This discrete field will serve as the analog of the Minkowski scalar 
field (/)(x, C). Notice that laboratory clocks tick off conformal time, so the time-dependence of 
these modes is as it would be for modes in Minkowski spacetime, and the time coordinate is 
appropriately labeled by C- Because the coupling matrix A only approximates a simple "balls on 
a string" coupling (characteristic of a bosonic field), the field will have a nonstandard dispersion 



relation, given by the frequencies Up from Eq. (17), as a function of the normal- mode index p. This 



deviation from ideality means that the analogy with a massless scalar field is an imperfect one. 



Prefactor and phase conventions for §„ vary in the hterature. 



B. Laser-induced coupling of vibrational and electronic states 

A feature of this model is that the internal electronic state of each ion can be turned into a local 
detector for the displacement of that ion. To achieve this for a given ion, its internal state must 
become correlated with its linear displacement from equilibrium. This correlation is provided by 
an external laser, which is used to drive an electronic transition between two meta-stable electronic 
levels \g) and |e), separated in energy by Hua [51 [T0l[22]. 

The interaction between an external classical laser field and the mth ion is described, in the 
dipole and rotating- wave approximationFlby the interaction-picture Hamiltonian [SlIlOj 



rC"*)/ 



M(/-\p+i[kcoseqmiO-uii^(;] _ A''''''-) ( ^\p~i[k cos BqUQ-uJLC] 



where Oq is the Rabi frequency for the laser-atom interaction (typically, Qq ~ 100 kHz [8]), k is 
the magnitude of the laser's wave vector k, which makes an angle 9 with the trap axis, (?m(C) is the 
interaction-picture position operator for the mth ion, and a^ {Q = e '^^'^a^ are its interaction- 
picture electronic raising and lowering operators, where a\^ = \e)^jj^g\ and a_ = \g)mmi^\- 
We may additionally absorb the laser's time-dependent phase e^^'^'^^ into the interaction-picture 
operators ^'"^(C) := e^'^^< a^JT^C) = e±^^^4"'^ where 

A = WA - WL (21) 

is the detuning of the laser below the atomic transition. This allows us to write the interaction 
Hamiltonian as 



r(™)/ 



M,'A\p+Jfccoseq„(C) _ Jrn),^. ikcoseq,n(0 



HPiC) = -ihno <j^-^(c)e+*'=-^^'?'"K)_^v_-^(^)e-*fecosyg„(g ^ (22) 

which is the common form in the literature [El dOj • 

The size of the rms fluctuation in qm as compared to the wavelength of the laser is measured 
by the Lamb-Dicke parameter 

'^■=y^M-u ^i^^^ose), (23) 

where ^Jh/2Ml' is the ground-state rms fluctuation of the center-of-mass mode of the ions. This 

quantity is representative of the overall rms fluctuations of the TTith ion, denoted Axrms, since the 
frequencies Vp for all higher normal modes of the ions remain within an order of magnitude of u for 
small numbers of ions (up to A^ ~ 10), realistic for current experiments [lOj . Typical values of the 
Lamb-Dicke parameter are 77 ^ 0.01 to 0.1 [23]. When rj^/fi + 1 <C 1, the so-called "Lamb-Dicke 
limit" (n is the average phonon number), the ion is well localized with respect to the wavelength 



of the laser, and we can expand the exponentials in Eq. (22) to first order in ry, which is equivalent 
to first order in k cos 6 qm{C,) , giving 

Hf'^ (C) ^ hno^ir^ (C) + hnok cos e g^ (C)4™^ (C) , (24) 

^^ — V — ■' ^ V ' 

carrier sideband 



There are actually two different rotating-wave approximations that will be mentioned in this paper. The first, 
discussed only here, is between the laser frequency and the electronic levels of the ion. This one is always valid 
since the laser frequency uji, ~ lja 3> v- References [S] [TD] discuss this in more detail. The second rotating-wave 
approximation — which is what will be meant by all further instances of this term — is between the detuning in 



Eq. (21 1 and the normal-mode frequencies. The validity of this approximation depends on these two frequencies, 



the coupling strength, and the time of the interaction. See Refs. [231 124) and citations therein for more details. 



where ^^'"^(C) = e+^^^af ^ + e-^^^ai'"^ and 4"^(C) = -ie+'^^a^^^ + ie-'^^a^^^. The first term 
corresponds to excitation of the transition directly by the laser, while the second couples the atomic 
transition to vibrational motion. 

With the laser sufficiently detuned away from resonance (|A| > z^) and having a sufficiently 
narrow bandwidth (interaction time longer than i^~^), we can neglect the carrier term and keep 
only the sideband term: 

^ uir) 

This Hamiltonian can be used to couple the atomic transition to the vibrational motion, with the 
detuning A serving as the effective resonant frequency for the detector. When A = fp (called 
the first red sideband), this interaction models a detector that becomes excited when a phonon 
is absorbed — an ordinary phonon detector. When A = —i/p (called the first blue sideband), it 
becomes excited when a phonon is emitted — a rather unorthodox detector indeed. After these 
interactions have completed, readout (projection) of the electronic state is accomplished using a 
bright laser that strongly couples the excited state to an auxiliary level, causing excited ions to 
fluoresce, while ions in the ground state remain dark |10j . 

As long as the coupling strength Qqi] is weak enough, no more than a single \g) — )• |e) transition 
will be excited for a given ion in the time during which the laser coupling is active. This allows 
us to use perturbation theory to calculate the detector response function. Strictly speaking, an 
infinite interaction time is needed for the simplified detector description given above to be exact; for 
finite times, the nonresonant terms (such as aJo"+) can be nonnegligible, resulting in red-sideband 
detections even when the ions are in the collective ground state. But this is what we want: this 
interaction has the form of a De Witt monopole coupling [21], and keeping these nonresonant terms 
is important for modeling the response of the modulated detector used to simulate the effects of 
expansion, the details of which we now discuss. 

C. Time-dependent laser coupling 

The usual construction, described above, produces a De Witt monopole detector for the ions' 
vibrational quanta. For our purposes, though, we need to introduce time dependence into our 
detector. What we need, specifically, is a detector whose Hamiltonian is scaled by a[t(C)]. But we 
can't simply apply this directly to the laser interaction Hamiltonian because the effective energy 
gap of the "detector" is hA, which is the difference between the atomic transition energy and the 
laser energy. Therefore, what we need is to apply the following modulation: 

A^A(C):=a[t(C)]A, (26) 

which, if we wish to vary only the laser frequency, can be achieved by 

UL ^ ujl{C) ■■= iOA - a[t{C)]A . (27) 

Alternatively, one could choose to Stark shift the atomic frequency such that cja — ^ '^a(C) ■= 
(x^L + a[t((^)]A. This may be more experimentally viable for some applications, but we will assume 
from now on that all modulation is done within the laser. 

The only other change necessary to simulate an FLRW metric in the detector picture is a new 
window function. Window functions correspond to modulation of the Rabi frequency fio of the 
atom, which is proportional to the amplitude of the interaction laser |10j . Since a window function 



10 



was needed anyway for initiating and ceasing the laser interaction, the requirement of a new window 
function presents no additional fundamental challenge. These modifications result in the following 
interaction Hamiltonian: 



U' 



Hf\0 ^^^^^^^ h^^r^FiO Y, ^(e-'^'-^a, + e'"^^ aW'^"-^ af^ + e-^^*(C)^(-)) . (23) 



(A^O) 



=1 /^r 



One other relevant thing to note is that the Lamb-Dicke parameter 77 depends on the laser wavenum- 
ber fc, which is proportional to uji^{C,). Thus, modulating the laser frequency will also modulate the 
coupling strength automatically through a time-dependent r/(C): 

7] k Ui^ WL 

but since wl ~ to a ~ 1 GHz, while A ~ z^ ~ 1 MHz, the fractional change in the coupling strength 
due to the frequency shift will be ~ 0.1% unless the scale factor gets to be of the inverse of that 
same order. This is unlikely, however, because such a large modulation (beyond a few tens of 
MHz) is both experimentally challenging [26l |27| and unnecessary because the detector would be 
so far detuned from the normal mode frequencies that it would (at best) detect nothing and (at 
worst) begin to excite extraneous atomic transitions. Therefore, while the time-dependence of the 
Lamb-Dicke parameter must be considered in actual experiments, it can be easily compensated 
by suitably modulating the laser power. We also note that a similar argument may apply to the 
Rabi frequency fio, depending on whether the atomic transition being excited is dipole-allowed or 
dipole- forbidden; in the case of the latter, Oq also varies with the atomic transition frequency [10] . 



which is of the same order as the laser frequency, and thus the calculation in Eq. (29) applies 



to ^{C), as well. Regardless of the experimental methods needed to achieve it, the new interaction 



Hamiltonian remains as in Eq. (28). 

At this point, we should take some time to talk about the approximations that could break 
down because of this additional modulation. One of the biggest concerns is avoiding unwanted 
atomic transitions. If the the bandwidth of the modulation, which is governed by the inverse of 
the detection time (width of the window function), is too large, this broadband response could 
overlap with unwanted transition energies. In addition, the range over which the frequency varies 
could cause the same problem even if the modulation is slow and over a long period of time. 
This issue strongly depends on the spectrum of the actual atoms used. Typically [8], however, 
v ^ 1 MHz, while uia ~ 1 GHz, which should provide quite a large bandwidth over which the laser 
frequency may be swept. Finally, the rate of modulation affects the feasibility of the scheme, as 
well. Modulations up to ~ MHz can be achieved with an acousto-optic modulator, while faster 
rates (up to ~ GHz) can be achieved using an electro-optic modulator [261 EZJ- We now describe 
several possible applications of our FLRW analog system. 

IV. EXPERIMENTAL APPLICATIONS 

A. Two distinct mechanisms by which expansion creates particles 

Birrell and Davies [18j go to great lengths to distinguish between two effects that might be 
called "particle creation" due to expansion. The first is the nonzero detection rate of an ordinary 
detector in what should be considered an "empty" universe when the detector is on during the 
expansion. If the expansion eventually stops, and the universe is asymptotically Minkowski, this 
effect is no longer seen. This does not necessarily mean that particles are not detected in that case, 
just that this particular mechanism for particle production only applies during expansion. 
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There is another type of "particle creation" that might be observed after a finite period of 
expansion has passed. Imagine a universe starting asymptotically in the Minkowski vacuum for 
very early times. Now imagine that expansion begins and ends smoothly, leaving the system in a 
nonequivalent vacuum but in a spacetime which is also asymptotically Minkowski for late times. 
Particle creation occurs in this case because of the Bogoliubov transformation mixing annihilation 
and creation operators of the original vacuum into operators of a single type in the new vacuum 
(oafter = ctObefore + /^Obefore)' ^^'^^ ensuring that there will be real particles flying around after the 
expansion. 

As we have formulated it, the detector picture can model the flrst case (detection during ex- 
pansion) but not the second (detection after expansion). The reason is the need for conformal 
invariance|j We have explicitly required that the vacuum we model in the detector picture be the 
conformal vacuum, which is (by definition) the ordinary Minkowski vacuum in conformal time. 
If we were to model a finite-time expansion in the detector picture, particles would be detected 
during the expansion, but none would be detected afterwardjj 

An ion-trap analogy has been proposed for each of these effects. Reference [TB] attempted 
to model the thermal response of a detector in de Sitter spacetime (exponential expansion) by 
modulating the trap potential and thereby modifying the phonons' normal modes. Since the 
proposal involved simulating the expansion by modifying the dynamics of the ions, the proposal uses 
the field picture. There were several problems with the proposal, including incorrectly predicted 
dynamics, calculational mistakes, and ions that would fiy out of the trap in a realistic experiment. 
These problems were identified and discussed in Ref . [T7] (see endnote [5] therein) , which itself was 
an alternative proposal for witnessing particle creation in an ion trap — again by chirping the trap 
frequency — but this time the proposal was for the second type of particle creation in which real 
particles are created by the simulated expansion and detected after it has ceased. This proposal 
does not suffer from the conceptual, calculational, and experimental problems of Ref. [16], but it 
also does not model the nonzero response of a detector during expansion. 

With the power of the detector picture, however, we can model a detector that is on during 
expansion — the original goal of Ref. [16] — but with a completely static trap, instead embodying the 
entire effect of the expansion within the modulation of the laser used to couple the ions' vibrational 
motion to their electronic states. In fact, we are not limited to de Sitter spacetime in this picture; 
we can model any scale factor as long as the required modulation is feasible experimentally. 

B. Single-detector response 

The detector response function Am for the mth ion, to lowest nontrivial order in the Dyson 
series expansion, is given by 

dC / < F(C)F(COe-^^W'^)-*(^')l {cl^UOMC')) . (30) 

-oo J —oo 

This is the probability that an ion will get fiipped to its excited electronic state by interaction with 
the field. The interaction is assumed to be weak enough to only excite one transition from \g) to |e) 
during the detection, thus allowing the use of time-dependent perturbation theory [28]. We have 
extended the limits to ±00, with finite limits able to be enforced using the window function F{Q, 



Our work does not rule out a detector picture for cases that do not admit conformal invariance. Our present 
formulation of it, however, requires this invariance. 

To see this, just imagine a detector whose resonant frequency is first fixed at one value, then smoothly shifted 
another value and then fixed at that value, but only after this process is the detector activated. For long enough 
detection times (much longer than the inverse detector frequency), the rotating- wave approximation may be made, 
and no particles will be detected [231124) . 
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as required. Also notice that the last term (in angled brackets) is just the ion analog of D^{C, ('). 



Since the field is the ions' collective ground state, we can plug in Eq. (19) and evaluate 



{cpm{C)<Pm{C') 



N r,{p)n2 

p=l 



(31) 



givuig 



N r,(p)i2 



p=i V^^p J- 



dCF{C)e 



'*At(C)g-ifpC 



(32) 



Notice that we can change integration variables back to t to obtain 






dtf{t)e 



-iAt—iUp({t) 



(33) 



This formula assumes n = 2, which is the case for a (l+l)-dimensional spacetime — this being the 
most natural spacetime to simulate in a linear ion trap, which has one effective spatial dimension. 
We note, however, that one-dimensional effects in spacetimes equipped with more than one spatial 
dimension can be simulated, as well, at the expense of additional powers of a{t) being included 



[see Eq. (14)]. The second form, Eq. (33), may be interpreted as the field-picture calculation of 



the same quantity as in the detector-picture formula, Eq. (32), since in this integral the detector 



time dependence appears ordinary, being just e~ , while the normal modes now experience time- 
dependent modulation, as e"*'^''^'*-' . Also notice that in the trivial case of no simulated expansion, 
C = t, and this quantity gives zero in the limit where /(t) — )• 1 (unless A < 0, indicating a blue- 
sideband detector), as a result of the rotating- wave approximation being applied to a detector 
in the vacuum. We can, however, still achieve a nonzero response in this case for short enough 
detection times [231 121] • 



C. De Sitter spacetime and the Gibbons-Hawking Effect 

One of the most celebrated features of de Sitter spacetime is the thermal response of a local 
detector to the conformal vacuum with a temperature proportional to the expansion rate: 



Tgh 



K 

2^ 



(34) 



in natural units. This is the Gibbons-Hawking effect [5]. It means that a single inertial detector in 
the conformal vacuum of de Sitter spacetime will respond as if it were bathed in thermal radiation 
at a temperature Tqh in Minkowski spacetime. An important physical difference between these 
two cases is that in the case of Gibbons-Hawking radiation, every inertial detector believes itself 
to be at rest with respect to the thermal bath it responds to — i.e., no velocity-dependent Doppler 
shift is observed [5] — while a thermal state in Minkowski spacetime has a preferred rest frame, 
and detectors in relative motion with respect to this frame (even if still inertial) will perceive a 
Doppler-shifted spectrum [18J. 

It is this thermal response that the authors of Ref. |16j wanted to simulate by chirping the trap 
potential and modifying the phonons' normal modes. Since the proposal involved simulating the 
expansion by modifying the dynamics of the ions, it uses the field picture. Later work corrected 
this proposal for a single ion [29^, which again involved exponentially modulating the trap fre- 
quency. The results were different from those calculated (incorrectly) in Ref. [16] but still showed 
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an effectively thermal response with a perceived thermal state with temperature TghJj Here we 
will instead propose to model the Gibbons-Hawking effect in the detector picture. 

De Sitter spacetime is an FLRW universe with an exponentially increasing scale factor in t: 

a{t) = e^' ^^ a[tiC)] = -^, (35) 

with transformations to and from conformal time given, respectively, by 

CW = -^ ^ .(0 = -!^. (36) 

which means that small intervals in each variable are given by 

dC = e-''^dt ^^ (it = -^, (37) 

where we have chosen i = to correspond to C = — 1/k. (In addition, t — )• oo corresponds 
to C — >• 0~, and t — )• — oo corresponds to C — )• — oo.) As in all simulations we are proposing, the 
clock on the wall in the laboratory ticks off seconds of conformal time C,, while the simulated proper 



time r of the detector is equal to the simulated cosmic time t and given by Eq. (36) in terms of C- 
Ideally we would like our detector to have an infinite detection time in t to get a perfectly thermal 
spectrum. But this is not physical, so we will simulate it with a detection time that is sufficiently 
long. Following the procedure proposed in Ref. [29j (but modulating the detector instead of the 
trap frequency), we will model the physical detector as having a flat window function over a flnite 
detection time from to T: 

f{t)^e{t)-e{t-T), (38) 

where Q{u) = 1 for positive u and for negative u. The equality is not strict because we want 
smooth turn-on and turn-off for the detector. The details of these functions is not important 
because the approximations we use allow us to model this detector as being on for a very long 
timeU] 

zv > K and i/e~'*^ < k . (39) 



A look at Eq. (33) reveals that in the field picture, the time-dependence of the field modes in 
an FLRW universe satisfies e~^^p^^^' . These limits, respectively, correspond to an expansion that 
begins in a regime where the rotating-wave approximation is valid (at t = 0) but extends into a 
regime where the field modes have effectively stopped oscillating (at t = T), meaning that particles 
will be detected j23[ I24j . The first approximation lets the lower limit of the window function be 
extended from i = to i — )• — cxd because the detector will see nothing during t < anyway. 
The second allows the upper limit to be extended from t = T to t — )• oo because the detector, 
which is assumed to be tuned to a frequency A ~ j^, is far off resonance from the mornal-mode 
frequencies Vp, which are all within an order of magnitude of v |10j . 

In the detector picture, this corresponds to sweeping the trap frequency from A (the actual 
detector frequency we want to simulate) at C = — 1/k to Ae at C = —e~/K, with intermediate 
values determined by 

A(C) = -|. (40) 



^ This is measured by taking the ratio of the detector response functions at the first red and the first blue sideband, 

giving g-^'^^/"-^ which is the desired ratio. 
'' We use the approximations from Ref. [29] rather than those of Ref. [16] because those of the latter (which correspond 

to rapid expansion of the trap initially) do not actually allow the limits on the window function to be extended as 

the authors claimed. 
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This is a polynomial modulation in (, while the proposals in Refs. [29] and |16| require exponential 
modulation of the trap frequency. The window function in conformal time (assuming n = 2) is 

^^(0 = -^. (41) 



Plugging these into the integral of Eq. (32) gives 



fO /-O 



/ dCF{C)e-'^'^^^e-'^^^ = f dCi-KCY^-^e-'^^^. (42) 

J —oo J —oo 

Substituting u = —k(, du = —ndC,, (3 = A/k, and a = Vp/n, along with appropriate convergence 
factors (a — )• a + ie, /3 — >• /3 — ze, with e — >• 0"'"), gives 

ird„u-''-'.- = -IM_ = !<!«,-«.. (43) 



Plugging in for the substituted variables and plugging the final result back into Eq. (32), along 
with the identity |r(z/3)| = 7r/[/3sinh(7r/3)], gives 

(n„,)^ 2. ^ [bi'f 

p=l V f 

which is a Planck spectrum at the Gibbons-Hawking temperature. This formula may be found in 
Ref. [16j, but it does not correspond to the procedure that the authors proposed. It is also different 
from the result obtained in Ref. [29] for the case of a single ion. This is expected, though, since 
modulation of the trap frequency creates single-mode squeezing that has no analog in de Sitter 
inflation. 

There are several experimental hurdles to be overcome. First, we want n to be low so that we 



can use the approximations in Eq. (39), but we don't want it to be so low that detection events 
are prohibitively rare (< 1%). Second, the laser modulation rates must be within feasible limits. 
Third, the electronic structure of the ions must be such that extraneous transitions are not excited 
while modulating the laser. We can address the first two problems briefly here while noting that 
any real application will require a more detailed analysis involving the atomic spectra. 



Combining the approximations from Eq. (39) gives e"**^ <C 1, which means that e"*^^ ~ 0.05 or 
less. This means that kC,{T) ~ 0.05kC{0), and so CiT) — C(0) ~ k~^, meaning the detector has to 
be on for roughly the inverse of the expansion rate. This makes the bandwidth of the modulated 
laser spectrum ~ k <C i^, so the normal mode frequencies are still well resolved, and we don't 
have to worry about broadband excitation of extraneous atomic transitions. The detector energy 
gap A(C) (and thus the laser frequency cjl) must be modulated over a range {e'^'^ — 1)A ~ 20 MHz 
for A ~ 1/ ~ 1 MHz [8]. The laser power similarly needs to increase by a factor of ~ 20 over the 
course of the simulation|£] 

While all of the above is achievable with current technology [261 [2Z! > the detection probabilities 
are prohibitively low due to the exponential scaling with A/k, which we have required to be large 



in Eq. (39) (since A ~ i/). Nevertheless, having established the theoretical connection to Gibbons- 
Hawking radiation for always-on detectors in our model, higher detection probabilities can be 
achieved using window functions of shorter duration and/or by using larger values of n. The 
resulting response function will no longer be a Planck spectrum but can be calculated numerically. 

* Simulating a higher-dimensional spacetime will change this requirement; see Eq. (fT4|. 
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This situation still models the analogous de Sitter detector — ^just one with a finite-duration window 
function. 

There are two goals, then, for achieving an experimentally accessible simulation. The first is to 
make the absolute excitation probability (for the red sideband) high enough to be detectable; the 
second is to have some sort of signature of thermality. To boost the detection rate, we will have to 



violate the approximations from Eqs. (39). This means the limits of the integral in Eq. (42) will 
be finite: 



C/ 



dCFiOe 



-i^t{C) ^-iUpC 



dnu^^-^e*"" 



T{i(3, —iauf) — r(i/3, 



-iauil^_^f3/2 



Ka'i 



(45) 



where u 



(iJ) 



-<{ij) 



e'^h^'f), 



and r(2;, b) 



J^dxx' 



^e ^ is the upper incomplete gamma 



function, which is a built-in special function in many computer algebra packages. We define the 
time interval in terms of the simulated initial and final cosmic times — ti and t/, respectively — in 
order to emphasize that this deviation from a perfect Planck spectrum is not an experimental 
artifact but results rather from simulating a detector that is on only for a finite duration (specif- 
ically, from ti to tf) in de Sitter spacetime. Defining the regularized upper incomplete gamma 
function Q{z,b) := T{z,b)/T{z), we can see that Eq. (44) gets modified to 



Am{A;ti,tf) 



(nov? 



2tt 



kA e^-"^l>^ - 1 



RiA;ti,tf), 



where 



N r;,(p)l2 
m 



R{A;t,,t,):=J2^^ 



Q 



A 



-i^e 

K 



-Ktl 



Q 



A 



-i-^e 

K 



(46) 



(47) 



which can be evaluated numerically. Extending the limits to ±00 gives Am{A) = Am{A; —00,00) 



in agreement with Eq. (44) 



In the case of effectively infinite-duration detectors, we can take the ratio of the red- and 
blue-sideband detection probabilities and find the signature of thermality: Am{A)/Am{—A) = 
g-27rA/K ^20j. For finite detection times that are shorter than required for the detector to thermalize 
with the perceived Gibbons-Hawking bath — which is the regime we are required to work in — this 



relation no longer holds. However, a new ratio can be calculated from Eq. (46): 



Am{A;ti,tf) 

Am{-A;ti,tf) 



-2nA/K 



R{A;ti,tf) 
R{-A;ti,tf) ' 



(48) 



which depends on when the detector is switched on and off but once again agrees with the fully 
thermalized result in the limit of long detection time. 

D. Multiple detectors 



The ion-trap analogy using the detector picture carries over to the case of many simultaneous 
detectors without modification. This is because all terms in the Dyson series are integrals involving 
the interaction Hamiltonian, and such integrals can always be recast from cosmic time to conformal 
time, resulting in a new interpretation as the response of Minkowski-spacetime detectors with time- 
dependent frequencies, as long as the window function is also appropriately modified. This leads 
us to suggest the possibility of studying the two-point correlations of scalar-field modes in de Sitter 
spacetime by measuring the correlation function for simultaneous electronic excitations in multiple 
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ions |28| under the detector modulation prescription proposed here. Such correlations are believed, 
in inflationary cosmology, to be responsible for structure formation in the early universe |13t [T^ . 
Two-component Bose-Einstein condensates have served in recent work [31j as analog systems for 
studying inflationary cosmology. Because of the restriction to conformally invariant settings in the 
current proposal, our goals are comparatively modest. By measuring correlation functions in the 
detection probability for excitations in multiple ions [28j under a detector modulation corresponding 
to various scale factors, we can analyze the spectrum of fluctuations that should result from the 
corresponding cosmological predictions [13l [U] . 

Independently of the cosmological implications, Ref. [32] proposes a method of swapping ground- 
state entanglement of a chain of ions to the electronic states of the ions at the ends of the trap 
using sonically local interactions (i.e., detection events that happen faster than the speed of sound 
can propagate across the trap). This is the ion-trap analog experiment of the proposal in Ref. |33| 
to use the Minkowski vacuum of a quantum field to generate entanglement between distant local 
detectors through entirely local interactions with the field. Using the detector picture, we can 
extend the proposal of Ref. [32] — as well as any other proposal that seeks to model only the 
Minkowski vacuum — to an equivalent proposal that instead models any FLRW universe in the 
conformal vacuum. This is particularly exciting because it brings into the realm of possibility 
performing the analog of the experiment proposed in Ref. |l5j to distinguish heating of a Minkowski 
field from de Sitter expansion through comparing whether entanglement is generated between two 
local detectors. The detector interactions used in Ref. [15] are impractical both in the actual 
vacuum of outer space and in an ion trap because the detection probabilities would be far too low. 
However, in an attempt to implement this analog experiment, the nonstandard detection techniques 
of Ref. [32J, which help boost the entanglement significantly over that achieved by the ordinary 
monopole detectors, may be applied to the two cases of (a) modulated energy levels and window 
functions, corresponding to de Sitter expansion, and (b) a heated field with no detector modulation. 
We leave the details of this proposal to future work, but we include it in the discussion to illustrate 
the wide-ranging applicability of the detector picture in ion-trap simulations of quantum fields. 



V. CONCLUSION 

We have demonstrated an alternative "picture" for analog QFT in curved spacetime, which 
encodes the effect of curvature in the detector rather than in the motion of the analog quantum 
field, and we have used this technique to propose new and more accessible analog experiments for 
curved-spacetime QFT effects. These include the Gibbons-Hawking effect for a single detector in 
de Sitter spacetime, plus the possibility of modeling structure formation in infiationary cosmology 
and the entanglement signature of de Sitter spacetime. 

A question which may be raised in the mind of the reader is the degree to which the ion 
trap/curved spacetime analogy continues to hold when time has been re-defined in such a way 
to make the trap resemble Minkowski spacetime. As one example, ordinary Minkowski spacetime 
has no cosmological horizon — one of the most conceptually important ingredients of, for example, 
de Sitter spacetime. Careful consideration of points like this show that the analogy still holds and 
that the horizon is actually still present. The horizon in the detector picture is related to the fact 
that evolution from — oo to -|-oo in de Sitter cosmic time corresponds to an evolution from — oo 
to in the analog Minkowski time variable C. The cosmological horizon is thus re-interpreted as 
the backward light cone extending from the detector location at (^ = in Minkowski spacetime, 
which equally well shields the detector from some of the field modes which are necessary to define 
the vacuum state. 

In a more general sense, however, all analog experiments are, at their root, quantum simulations. 
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In a simulated system, there is no reason to expect that simulated time must evolve in the same 
way as ordinary time. Anyone who has done numerical simulations of time-dependent systems can 
attest to this. Going between the two pictures can be thought of as running exactly the same 
simulation program on two different machines — while the simulation may be run at different rates 
(depending on the machine), the end result is precisely the same output, and the degree to which 
the output is analogous to something in the "real world" is identical. 
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